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THCUNICAL MEMORANDUM 


STATIC AND DYNAMIC DEFLECTION STUDIES OF THE SRM AFT 

CASE/NOZZLE JOINT 


I. INTRODUCTION 


A. Background 

Recently, there has been an increasing interest in efficiently developing large math model 
eonl.gurat.ons for. dynamic and static analyses. The prediction of the behavior of a structure is 
genet ally based upon the results ol the analysis of a mathematical model of the structure The 
accuracy of this prediction depends on how well the model approximates the structural characteris- 
es. Hoi this reason it is important to know the limitations of modeling procedures, the software, 
and (he model developed to represent the structure. 

I he rapid development of computers has completely transformed research and engineering 
practices m every technical field, particularly in the field of structural dynamic analysis. Conven- - 
(tonal computers have rapidly branched out into supercomputers, minicomputers, and graphics 
computers and workstations. Personal computers today are as common as calculators were in the 
U7() s and sltde rules '» *he l%() 's before. Microcomputers are being developed with 
unprecedented memory, speed, and graphics capabilities. Most of the math modeling experience 

has been accumulated m the last two decades primarily because of the development of and increase 
m the use ol computers. 

Math model characteristics should be similar to that of the physical structure. Models are 
used to p an. design, and study the physical structure. In most cases, modeling and analysis of the 
model reduces the cost, risk, and the amount of time required to accomplish these tasks. Frequent- 
ly. math modeling ol the structure is the only way to do this because the hardware may be 
unusable or unavailable for testing or for other purposes. 

In order to resolve problems using math models, it is necessary to understand the system (or 
s ructure) and the problems relating to that system. Math models should be developed, from the 
stait, with a specific problem or set of problems in mind to facilitate best use of the model. The 
software used to translate the model into a form usable by these newer generation computers must 

, bc “ su ' lab,c veh,cle to P erform this ‘ask. The results (output) should also be able to he 
analyzed and correctly interpreted such that engineering decisions to resolve the problem can be 
made. I he success of a modeler depends on how well he/shc can define significant elements of the 
pioblem, establish the relationship between each part, and obtain meaningful and accurate inlorma- 

; n '" ,n !‘ e mt ; del - , Th ; S re u p ° rt prcsents a technical slud y of this 'W for the solid rocket motor 
(SRM) alt domc/nozzlc lor the Space Shuttle vehicle. 


B. Purpose 

CDM ? evem ! typcs of structural dynamic and static response analyses were performed on the 
a,t case /nozzle joint, with ali-thc results from these studies being presented in this report, 
te main purpose of this study was to determine -gap deflections at the joint caused by the 
chamber pressure, nozzle blowout forces, actuator gimbal loads, and- vehicle low-frequency vibra- 

lon-accelerations. The computed data is compared with test data to verify the model of the SKM 
structure » 


The location and configuratiou of the joint can be seen on Figure I with the enlargement on 
igure 2 showing a close-up- of the joint with the new design using the- radial bolts, (The old 
design did not include the radial bolts.) The aft dome, fixed housing, and nozzle were modeled 
using MaeNeal Schendler NASTRAN (MSC/NASTRAN). Figure I shows a cross-section of these 
parts and their relative positions. The model was fixed in the X, Y,. and Z directions at the tangent 
point, the tangent point being located a short distance below the aftmost factory joint on the aft 
segment of the SRM as shown in Figure I. 


TANGENT 
POINT 

(SRB STATION 
X- 1823.95) 
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Figure I. SRIl aft dome/nozzle joint. 
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Figure 2. Current joint design. 


II. MODEL DEVELOPMENT 


The analysis of the SRM aft dome/nozzle joint structure involved developing a large finite 
element model of the aft dome, the fixed housing, and the nozzle. The model was developed using 
an Intergraph Graphics Design Station (1GDS) workstation for input to MSC/NASTRAN software 
on the IBM/Cray computer. The structural dimensions and details were obtained from Morton 
ThiokoL drawings noted in References ! through^ 

An isometric view of the model is shown on Figure 3, which started with approximately 
32,400 degrees-of- freedom (DOF) and then was reduced to approximately 22,050 DOF. The nozzle 
is modeled with beam elements and includes two rod elements for. the nozzle actuators. Plate, bar, 
rod, and solid elements are used to model the complete joint structure. The aft skirt was not 
included in the math model because it is much sliffer than the aft dome, to which it is connected, 
so that the analysis would provide conservative results. The math model was developed in the full 
360-degrce configuration to allow for unsymmetrical- loading and checking results in an 
unsynimctricnl fashion to be sure the worst case was enveloped. 

The model was developed with more fidelity (finer grid mesh) at the fixed housing/aft dome 
interface (Fig. 2) with 2,100 grid*; between these parts. At the SRM case tangent point, the 
boundary/base of the model, the model has the coarsest grid spacing with 50 stations around the 
circle. This spacing compares to 200 stations around the circle at the aft dome/fixed housing 
interface. The fixed housing and aft dome arc held together with both vertical (axial) and radial 
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Figure 3. SRB aft dome/nozzle joint dynamic model. 

bolts. These bolts are modeled with three grids and two bars, one bar which has the bolt preload 
applied to it. This preload helps lower the gap opening when the loads, such as the SRM internal 
pressure, are applied to the structure. The effects of the bolt preloads are included in the analyses 
results and will be evident in the data presented. The 50 boundary support points were all 
constrained in the three translational directions. The low frequency vibrational motion at these 50 
points were available only in the three translation directions, therefore the rotational DOF's were 
removed from these boundary points. 


The flex bearing assembly in. the center area of the model and the vertical bolt area of the 
model (Figs. I and 2) were made with slightly less fidelity than the radial bolt area of the model 
o minimize the number of DOF and speed up the analysis time. For example, the stiffness of the 
Ilex bearing area is not directly included at the fixed housing/nozzle interface, but the nozzle part 
of the model was tuned to nozzle/llex bearing frequencies |5,6J to account for this llexibility in the 
model and match the structure's behavior. To help accomplish this step with the model, the rota- 
tional stiffnesses at the fixed housing/nozzle interlace was left out to represent the ilex bearing 
mom accurately. The nozzle is represented by a cantilevered beam with dynamic characteristics 
similar to a larger, more detailed and complete nozzle model, 

T he primary purpose of the model was to provide deflections along the vertical face at the 
locations of the primary O-ring and the radial bolt. The modeling procedures used in developing 
this model reflect this purpose. This report presents only the deflections along the vertical face 
between the alt dome and the nozzle fixed housing. Figure 4 shows a section of the aft dome/fixed 
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housing interlace with these locations shown as dots where gridpoints were on this vertical face. 

For example, a dot (also a gridpoint), located at station -50.587 inches is at the centerline of the 
radial bolts. 


III. THEORETICAL TREATMENT - ASSUMPTIONS AND EQUATIONS 


, . llK n 0q ‘ ,a,io,1s , r ° r this of study primarily involve the use of large mass and stiffness 
matrices. the _ Unite clement technique, using NASTRAN for example, is a way of constructing the 
mass : „,d sill ness malriees. These malriees are ihen used to develop a group of seeund-olr dlf 
kicnnal equations which are to be solved simultaneously. The equations also include damping 

sm!c,urc P !nT n,S abi,i| y 10 absorb en ergy. Many types of loads are applied to fhe 
stfUUurc and are included m these equations. The equations, including mass, stiffness, damping 
and applied forces, are as follows: 1 


|M| {X,} + |C| {XJ + |K) {X,} = {FJ 


where 
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[Ml — mass matrix 

[C| - damping matrix 

I K 1 — stifl'rtGss matrix 

{Fj} - forces acting at grid i on the model 

{Xj} - general vector representing the deflection at grid i on the model. 

The mass matrix is an array representing the weight of the structure distributed over the grid 
points of the model, while the stiffness matrix is an array representing the stiffness between the 
grid points. The force vector on the right side of equation (1) arc the forces acting at various grid 
points on the model. These forces include internal chamber pressure, blowout force, pressure forces 
and moments from nozzle gimbal, actuator gimbal forces, gravity, and qua, '-static load factors. 

The pressure forces are applied in several different ways. The chambe. pressure acts over the all 
dome, the fixed housing, the nozzle, and other areas of the SRM. However, the math model in 
this report includes only the aft dome, fixed housing, and nozzle. The primary area of interest for 
this report is the worst case gap deflection, at the aft dome to fixed housing interface around the 
entire 360 degrees of the model. The chamber pressure acting on the aft dome and fixed housing is 
applied as a force acting at each grid point. The pressures acting on the nozzle are integrated over 
the surface of the nozzle to establish forces and moments at the nozzle alt end ling face (Fig. 1). 
The gas dynamic equations used (from Reference 7) to establish the nozzle forces and moments 
were developed by the Morton Thiokol Space Division and are presented in the following for 
completeness. 

The axial com ponent (orbiter X-direction) of the force is: 


Fnuo = A n * P thn + T V ac (• - cos 0) + P Alllb (A,.;|, cos 0 - A|. ; i<) 


( 2 ) 


where 

Fnuo — nozzle axial blowout force - vertical component of gimbal nozzle force (parallel to 
motor axis), lb, 

A n - 1 ,026 in 2 - ratio of null position nozzle force to nozzle stagnation pressure, 

P c hn - nozzle stagnation pressure, psia, 

Tv Ac = vacuum thrust, lb, 

P Ami* = ambient pressure, psia, 

A|;|» = 'nil A. 822) 2 = 17,588 in 2 - nozzle exit plane area, 
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/'fin - ir {M.2.Z) = 4,832 iir area enclosed by afl end ring O-ring, 

0 ~ gimbal angle, degrees. 

Pic lateral component of the force (Ih) due to gimbal angle is: 

Fn,,„ = (Tvac — At-P P Amh) s i>l 0 

The force direction is opposite to that of the gimbal direction. 

The yaw/pitch moment (in-lb) at the aft end ring due to nozzle gimbal is: 

My/i'ii = Z|>|» Fnlu, 


where 

Z|.|» = 6.057 in = moment arm between nozzle pivot point and aft end ring face. 
The lateral component ol the transient force (lb) due to gimbaling rate is: 

FNUtTr = (4.521 x 10 - s ) Ln m 0 cos 0 

where_ 

L n = 175.757 in = nozzle nose to exit length, 

0 = 5 deg/s = maximum nozzle gimbaling rate, 
m = total mass flow rate,Jbm/s. 

The axial component of the transient force due to gimbaling rate is: 
pNnoTr — (4.521 x | O' 5 ) L N rh 0 sin 0 

The transient yaw/pitch moment at the aft end ring face due to gimbaling rate is: 


Mv/mv « (4.521 x 10 ^ Z N L N rii 0 



where 


Z N ~ 53.221 in moment arm betweeu-uoz/le midpoint and alt end ring lace. 

liquations (2) through (7) give the forces caused hy the gas dynamics involved while operat- 
ing the SRM. -These forces, as described, act at the nozzle alt end ring face with the associated 
moments, liquation (2) gives the nozzle blowout force which acts in the direction of the Space 
Shuttle Vehicle (SSV) X-axis on the nozzle and tends to push the nozzle downward, (i.e. , the SSV 
+ X direction), liquation (3) gives the lateral force acting on the nozzle caused hy thrust and 
nozzle gimbal. liquation (4) gives the moment acting on the all end ring face with the force Hni. UI 
also applied at the aft end ring face, liquations (5), (6). and (7) give forces cau the nozzle 
gimbal rate and the fact that there is a lluid mass moving through the nozzle. These forces and 
moments, in most cases, arc going to be small and will have little influence on the gap motion 
studied in this analysis, compared to the effect of -the internal chamber pressure and the radial bolt 
preloads which have more influence. 


Now that the various parts of equation (I) have been described, the solutions to this equa- 
tion and their uses will be discussed. The direct transient response analysis, the static analysis, and-- 
the base driving response analysis were alt solvcd-using different techniques which will be 
discussed in the following paragraphs. 


A. Direct Transient Response Analysis 

The direct transient response approach was used to calculate vibrational motion of the gup 
caused by the buildup of the SRM chamber pressure at ignition and gimbal forces (caused by the 
nozzle). This approach is referred to as a direct approach because equation (I) is used in 
NASTRAN and solved without calculating modes and frequencies of the model. In other words, 
the modal response method was not used. The equation was solved numerically in NASTRAN 
(i.e., not in closed form). The input forces build up from zero to their maximum values using the 
same time scale and shape of curve as the chamber pressure rise curve. The loads applied were the 
internal chamber pressure, the nozzle blowout force, and the side forces and moments from nozzle 
gimbal as shown in equations (2) through (7). The chamber pressure is distributed over the aft 
dome and lixed housing while the blowout force is acting on the top end of the nozzle and the 
gimbal side forces are applied at the nozzle alt end ring face. The results of this study are 
presented in Section V. 


B. Static Analysis 

The lorces used in the static analysis were the maximum forces from the direct transient 
response analysis while also including the nozzle actuator forces and the liftoff load factors. The 
static analysis was also conducted using NASTRAN. The stiffness matrix of the math model was 
used to obtain deflections and stresses. The mass matrix was also used in the analysis to apply the 
liltoll load factors at each grid point of the model to establish the acceleration induced loads. The 
nozzle actuator forces were applied to the nozzle at the locations where the bottom of the actuator 
ties into the nozzle. Ihese lorces were calculated using the maximum nozzle design bending 
moment of 4.2 1 ().()()() in-lbs (obtained from Reference H). 
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C. Base Driving Response Analysis 

I lie base driving response approach was used to determine the influence that vehicle 
dynamic motion (vibration) has on the gap del lections at the all dome to fixed housing interlace. 

I lie vehicle dynamic accelerations were taken from Rockwell loads studies performed on the entire 
shuttle vehicle (launch configuration). The accelerations were defined- in these studies at four 
positions around the case at SRR station 1 .606.72 in. 90 degrees apart in the SSV coordinate 
system Y-Z plane. The math model has its boundaries at station 1.K2.1.95 in. Therefore, one of the 
assumptions for this analysis was that the accelerations at the two SRB stations would be similar in 
boili magnitude and frequency content. These accelerations were available at these points around 
both the leti and righuSRU's. The worst case accelerations were obtained from these data in each 
o! the thiee directions and placed at all 50 boundary points around the model to perform the base 
drive analysis. The equation used to perform this operation is derived in the Appendix as equation 
(26) and is presented here as equation (8). 


[M«,l {q} 4 |2£u>| I M 0tJ | {q} 4- |co“| |M Ct ,| {q} = |(J)| T ({F 2 } - |M 2 ,| {X,} - [M-l_j£] {X,}) 


( 8 ) 

where 

|M C ,,| = |<|)| l | M _. j | |(|)| (generalized mass) 

|2£w||M cl| | = [(hi 1 |C— | |t|>| (generalized damping) 

U^illM Vl| | = | <1> | 1 |K-| |t|)| (generalized stiffness) 

and. 

|M-| = |M| 

|C-| ..= \C\ 

IK-1 = | K | 

are from equation (I) of this section. 

The other variables are 

{X|} - boundary base driving accelerations, 

{F : } - applied forcing functions = {()} for base driving, 

|(4| - rigid body transformation that changes boundary accelerations to body (model 
gridpoint accelerations). 
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,M -,I " P arliUm) ° r mass main* related u> the ground points; only-spears in the free-free 
partition ol the matrix, 

{q} - generalized modaLeoordi nates, 


N>l 


the modal deflections calculated Ironv equation (I). 


Thc How chart lor the computer program which runs this analysis is diagrammed in Figure 
5. The program resides on computer files and is run using FORTRAN programs. The chart. begins 
on the left with data on file created using NASTRAN and other FORTRAN programs. The chart 
ends on the right with response versus time of accelerations and deflections. It can he seen that- 1 ha 

inal deflections versus time in physical discrete coordinates are calculated from the following 
equation; 


{X 2 O)} = 101 {X} + |<|>j {q} 


where X is the motion at the support points of the hardware and q is the vibrational motion of the 

point m generalized coordinates. The product of |<J>| and {q} transforms the deflections into physical 
discrete coordinates represented by X 2 (t). 


IV. TYPES OF ANALYSES 


A. Modal Analysis 

The modal analysis was conducted on the aft case/fixed housing model with the nozzle 

beam model to obtain the modes, frequencies, and the mass and stiffness matrices for use in the 
base driving response analyses. 


A detai,ed madl m °del ol the ait case to fixed housing was developed to ensure a high 
degree ol lidelity at the joint, but the nozzle was represented by a beam model developed with a 
minimum number ol DOF. However, to ensure accurate results, the nozzle beam model was tuned 
o match results ol a very detailed modal analysis of a cantilevered nozzle, which included steps of 
tuning the model without the actuators, and then with them attached. This detailed analysis is 
published in References 5 and 6. 

The SRB mass property document was also used to obtain the correct mass for both the 
detailed part ol the model and the simplified nozzle part of the model as well as (he total mass for 
the combined model. I he masses Irom this analysis compared well with the masses found in the 
mass property report l he frequencies and mode shape descriptions of this model are shown in 
able I and the mode shapes are plotted in Figures 6 through 12. The lower frequency modes 
describe mo, ion ol the nozzle which is represented by beams and arc seen as lines coming out of 
the bottom center ol the model on the views shown of the modal plots. Modes 6 and 7 involve 
most ol the flexible motion of the entire aft dome and fixed housing. 
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Figure 5. Base driving analysis flow chart. 
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22.28 Hz 
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65.47 Hz 


185.5 Hz 
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of dome ° 
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Figure 6. Modal deformation: Mode I, frequency 2.22 Hz. 
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Figure 9. Modal deformation: Mode 4, frequency 22.28 Hz. 



Figure 10. Modal deformation: Mode 5, 


frequency 22.32 Hz. 



B. Static Analysis 


The sialic analysis involved applying (he maximum* lor every type of load considered. The 
loads applied were internal chamber pressure, nozzle blowout forces, actuator gimbal forces, verti- 
cal and radial bolt preloads, gravity magnified by a steady state acceleration factor, and vehicle 
low frequency vehicle dynamic accelerations applied at all the grid points in the model. This was 
done for each flight event, a total of five different subcases (flight times). Kach flight event 
involved a subset of these forces using the applicable values for the maximum in each subcase. 


The first event considered, liftoff, was the only case in which the vehicle low frequency 
vehicle dynamic accelerations (quasi-static) were applied. The other four events - roll maneuver, 
maximum dynamic pressure (max q), maximum acceleration (max g), and pre-staging (SRB separa- 
tion) - include the remainder of the loads considered. A table summarizing these values used in 
each of the five cases is shown in Table 2. 


in addition to all these cases, the model was run without the radial bolts, as in the old 
design, to compare the behavior of the no-bolt configuration to that of the new 7/8-in radial bolt 
fiight configuration. The force values used in the no-bolt analysis are also the ones shown in 
Table 2. 


C. Transient Response Analysis 

In the transient response analysis, the opening of the gap was computed as a response to the 
input pressure transient (Fig. 13) and other forces related to the pressure transient such as the 
blowout forces. This analysis was performed to observe the dynamic effects (overshoot, amplifica- 
tion. etc.) on the gap opening compared to the static value of the gap opening. For this reason, the 
slope of the chamber pressure rise curve is important. For example, if the chamber pressure curve 
rises to a peak from zero time to a time equal to one-lourth the period of the first natural fre- 
quency of the structure, then large vibrational responses in the structure can occur. 

An enveloping pressure curve was obtained from References 9 and 10 and used to force the 
math model. The curve is a smooth pressure rise curve with the dynamic behavior of the model 
mainly sensitive to the slope of the curve. The pressure rise curve builds up in 0.2 s. This would 
resonate with a frequency of approximately I Hz. 1/(4 x 0.2) = 1.25 Hz. The lower frequencies of 
the whole model with the nozzle are 14 Hz to 22 Hz and the natural frequency of the aft dome/ 
fixed housing is 65 Hz. It can be seen by comparing these frequencies to the forcing function fre- 
quency that there will probably be very little dynamic response due to the pressure transient curve. 

A check case was run with the slope of the pressure transient curve increased by a factor of 
10. The results show an increase in the amount of the dynamic response, but compared to the 
static pail of the deflection, the dynamic portion of the deflection is still very small. These 
response plots are shown in Section V. 
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TABLE 2. SUMMARY OF FLIGHT LOADS 


Internal pressure . 920 psj 

Steady slate 

acceleration (G's) 1 .4 


Blowout forces: 
Axial (lb) 

Side (lb) 
Moment (in-lb) 

Actuator loads: 

X 

Y 

Quasi-static load 
factors:* 

X 

Y 

Z 


1.13 + 6 
113.649 
24,653 


23,320 
-5 1 ,62 1 


+ 1 .0/-4.0 
± 1 .6 
± 6.8 


(Nozzle mass = 18,873 lb) 


:|i From Reference 8. 


Roll Maneuver 

Max Q 

Max G 

Pre-Stage 

930 psj 

788 psi 

673 psi 

90.7 psi 

1.4 

2.0 

3.0 

3.0 

1.14 + 6 

8.60 + 5 

7.07 + 5 

9.30 + 4 

285,261 

265.833 

171,362 

43,609 

25,616 

21,923 

19,157 

2,960 

23.320 

23.320 

23,320 

23,320 

-5L62I 

-51,621 

-51,621 

-51,621 
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Figure 13. Measured pressure on aft SRM dome. 










Dr Base Driving Response Analysis 

The base driving response analysis was conducted- using the model analysis data described 
above in Section IV. A. Vehicle acceleration time histories from coupled loads analysis were -used 
at the model attach points-to drive the model. Typical plots of these time histories at a grid are 
shown on Figures 14, 15, and 16. 

The How chart (Fig. 5) shows all the steps necessary to calculate loads and deflections from 
ase input accelerations. Steps I and 2 (step numbers are shown in the upper right corners of the 
blocks), show the accelerations and enforced displacements used to drive the model as input to 
start the analysis. Steps 3 through 14 show all the matrix manipulations performed using the mass 
and-stiffness matrices. Both the free-free matrices, M„ and K, r , were input to the. base driving 
program and disassembled into the ground motion partitions, M 2 | and K 2I , and the cantilevered 
partitions, M 22 and K 22 . It is also possible to obtain these cantilevered matrices directly from a 
modal analysis with the base of the model constrained. To make the partitioning step a single 
operation, the ground points in both matrices were repositioned to the first rows (and columns) 
while using NASTRAN. The K 22 partition of the stiffness matrix was inverted and combined with 
the other partition of the stiffness matrix to form a new matrix called the beta (fi) matrix. This 
new matrix is then combined with the M 22 partition and added to the M 2i partition of the mass 
matrix to form a coefficient matrix for the base input accelerations. This coefficient was developed 
in equation (8) of the previous theoretical section. 

The mode shapes, natural frequencies, normalized mass matrix, and the time vector (steps 
15 through 27), along with the newly formed forcing function matrix, are then fed into a modal 
response subroutine shown in Figure 5 as ZTR3. The modal responses from the ZTR3 routine are 
then used as input for another response subroutine (ZTRAE) where these responses are transformed 
back into physical domain accelerations, velocities, and deflections. These physical coordinates 
describe the model’s response to the original input forcing functions giving the gap deflections and 
loads as desired. The results of this analysis are shown in Section V. 


V. RESULTS 


A. Model Verification with Transient Pressure Test Article (TPTA) Test Results 

The TPTA test series were used to verify the aft case/nozzle joint model for these analyses 
by using the math model with slight changes in its configuration to model TPTA conditions 
accurately. Details ol the TPTA test setup were obtained from Reference II. The TPTA test used 
loads matching the flight values except for bolt preloads and nozzle loads. The test also had 
slightly dilferent boundary conditions. The actual test set-up had a plug fitted tightly into the fixed 
housing/nozzle interlace area with an O-ring sealing this joint for the TITA test. 


The main modification to the math model was removing the beam model of the nozzle atld 
fixing the boundary conditions at the fixed housing/plug interface in all translational directions 
except the X-dircction. This, most likely, makes the model stiller than the test article in this area. 
Other modifications included lowering the radial and vertical bolt preloads from the (light values to 
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Figure 16. Base driving forcing function Z-direction. 











the lest values. Typical prdo;uLv;diies for (he. flight configuration for the vertical and- radial holts 
are 141,000 and 45.000 Ih, respectively. The preloads- for the. TPTA- tests compare at 92,000 and 
25,000 lb for the vertical and radial holts. 


The deflection gauge locations used in the TP PA test were at the same circumferential loca- 
tions as the model where grid points and displacement information was retained. These locations 
were obtained from Morton Thiokol drawing No. 7U75234 |12|. 

The plot in Figure 17 shows the individual data points from the first TPTA test as stars and 
the analytical results from the model as circles connected- with a line. The two sets of data arc 
plotted together to show the similarity between, the test results and the analysis results. The analyti- 
cal results envelope the test data which shows that the model is conservative for studying the joint 
behavior. 


B. Analyses Results 

Results of the static and dynamic deflection analyses for the flight configuration are covered 
in this section. In general, the maximum envelope results of all cases are shown. Four locations 
with grid points were used in the model on the vertical face at the aft case/fixed housing interface 
shown in Figure 4. The results for the vertical face are summarized on Figures 18 through 22. The 
actual deformations calculated by the model along the vertical face are shown as deflected shapes 
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Figure 17. Static gap between fixed housing and aft dome - TPTA test. 
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Figure 18. Static gap between fixed housing and aft dome - flight niath model. 



Figure 19. Static deflected model compared with undcflected model. 
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1‘igure 20. Gap deflections with and without radial holts. 
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Figure 21, Transient analysis gap deflections. 
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Figure 22. Relative vibration gap motion at aft dome/fixed housing (X -direction). 

on Figure 18. Figure 18 shows the deflections at different flight events plotted to scale, while 
Figure 19 gives the cross-section deflections amplified by a large scale-factor so that the general 
shape of. the deflection pattern can be studied. 

Figure 20 shows two static deflection curves, in the radial direction, for the vertical part of 
the gap between the aft case dome and the nozzle fixed housing for 930 psi chamber pressure and 
4.8-degree nozzle gimbal angle. Curves are shown for the configuration without radial bolts (old 
design) and with the 7/8-iuch radial bolts (new design). Deflections at the centerline of the primary 
O-ring gap can be seen by following the horizontal axis to the markings on the graph and going up 
to both curves. The left end of each curve is the bottom of the vertical part of the gap between the 
alt case and the fixed housing and the right end is the top. As can be seen from the figure, the 
7/8-inch radial bolt reduces the gap opening by almost an order of magnitude with respect to the 
configuration without radial bolts. Figure 18 shows the results of the studies using the 7/8-inch 
radial bolts for different flight times and conditions (chamber pressure and nozzle gimbal angles). 
The top curve is for the vehicle roll maneuver case. This case was run using 930 psi as the motor 
interna! chamber pressure, a nozzle gimbal angle of 4.8 degrees, with blowout forces, actuator 
forces, steady-state vehicle accelerations, gimbal side force and moment. The next worst case is 
liftoff, with 920 psi chamber pressure, 2-degrees nozzle gimbal angle, liftoff vehicle accelerations, 
and the same type ol forces and moments as indicated above for the roll maneuver. The third 
curve is for the maximum dynamic pressure event with parameters of 788 psi chamber pressure and 
5-degrees gimbal angle. The last case shown is for the maximum acceleration case using 673 psi 
and 3. 7 -degrees gimbal angle. 

Transient response results during liftoff are shown as gap relative deflection versus time on 
Figure 21. The four curves are the four locations on the vertical face where the relative deflections 
were calculated as shown on Figure 4. These curves follow the shape of the pressure versus time 
with only a slight amount of vibration as shown on the curves in the period of time when the 
maximum pressure is reached. These curves are presented to show the small amount of vibration 
which occurs in the gap between the aft dome housing and nozzle fixed housing as a result of the 
applied transient forces. The maximum deflection at the top of the gap is 0.0063 inch as shown. 


Studies were also performed using liftoff vehicle dynamic motion base drive the nozzle/aft 
case joint structural model with acceleration time histories from STS coupled loads analysis |I3|, 
The results of this study show extremely small deflections at the nozzle/aft case joint gap as a 
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result of vehicle low frequency vibration, input. This can be seen on. Figure 22 and shows the rela- 
tive-deflection between hoth-sides at- the top of the gap. As can- be seen, th& relative deflections 
arc very small, compared to the other analyses (of order < l_x I O' 6 ). Other deflections down the 
lace of the vertical part of the gap will he smaller. 

These analyses have shown that-nozzlc-to-cu.se joint dynamic gap deflection motion caused 
by SRB internal chamber pressure buildup, nozzle gimbaling. and vehicle low frequency dynamic 
motion is very small compared to the static gap deflection. The maximum gap opening at the 
centerline ol the primary O-ring is approximately 0.0049 inch using the. 7/8-inch radial bolt con- 
figuration with a nominal preload of 45.000 lb. 


VI. CONCLUSIONS AND RECOMMENDATIONS 


A very large Unite element analytical math model of the SRM aft dome/nozzle joint for the 
Space Shuttle vehicle was developed to investigate gap deflections at the O-ring seal during simula- 
ted flight operation. All in-flight static and dynamic environments or forcing functions were applied 
to the analytical model to calculate the gap deflections. In addition, the finite element math model 
was modilied slightly to simulate an experimental program conducted to measure the gap deflec- 
tions. The analytical and experimental results closely match. This comparison provides confidence 
in analytically predicting the gap deflections during flight, and is considered a step forward in the 
use oi linite element modeling techniques lor large complex structures. It is recommended that this 


type study be conducted in the luture on all O-ring interfaces for high pressure solid rocket motors. 
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APPENDIX 


Response Equations for Base Motion Excitation 

By Wayne Holland, S&E*ASTN-ADL 
March I, 1971 


The equation of motion for a discretely modeled structure can be written in matrix form as 


[M] {X} + [C] {£} + [K] {X} = {F(t)} 


(1) 


where: 


IM], [Cl, [K] = mass, damping, and stiffness matrices, respectively 

{X}, {X}, {X} = acceleration, velocity, and displacement vectors of discrete coordinates. 


Consider that the system of equation (I) has n number of degrees of freedom. That is |M|, [CJ , 
and [K] are n x n matrices and, correspondingly, {X}, {£}, {X}, and {F} are n x I matrices 
(vectors). Additionally, consider that m number of the displacements are known or prescribed. The 
known or prescribed displacement functions are referred to as base motions. 

We can partition equation (I) to reflect the prescribed coordinates in the m x l vector, 
{X|}, and the remaining coordinates in the (n-m) x I vector, {X 2 }. With this partitioning scheme, 
equation (1) takes the form 


M,,jM 12 

1 

M 2I 'M 2 2 


X! 

--) + 

*2 


C,,jC, 2 ’ 

c 2 ,!c 2 2 



K,,}K 12 

j 

K 2 | iK 22 




X, 

X 2 


F,(t) 

F ? (t) 


( 2 ) 


Let {P} be the n x 1 vector of resultant external forces (includes applied and inertial 
forces). Then the relationship between resultant external forces and the displacements can be 
written in partitioned form as 



K||jK | 2 

- 1 


K 2 rK 22 

» m 

1 <N 

Li 


(3) 
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Now lei us define the {X 2 } displacements to be a superposition of two parts 

)X,( = {X^ + {X,} (4) 


where: 

{X 2 } = displacements that occur at the {X 2 } coordinates due to {X|} displacements, con- 
sidering the {X 2 } coordinates to be unloaded; i.e., with {P 2 } = 0 

{Xj} = displacements that occur at the {X 2 } coordinates considering the {X,} displacements 
constrained to be zero and with the {X 2 } coordinates loaded; i.e., with {P 2 } ^ 0. 

Now we can determine {X^} by applying equation (3) with {P 2 } = 0. 

Thus, 


p. 


~k,,[K I2 ' 

(X, 

0 


I . 

!2_j_ 

*T 
« ! 

I'xY 


It follows from equation (5) that. 


)x“} - fpi {X,} 


(6) 


where: 


ipi = - fK 22 r' (k 21 i 


(7) 


In view ol equations (6) and (7), we can now write equation (4) as 


{X 2 } = IP! {X,} + {X 2 } 


( 8 ) 


With equation (8) we can now write the transformation equation in the form 





» [07 jx , j 

pH (x 2 ) 


( 9 ) 


or more simply 


W = l7l {Y> 


where: 


( 10 ) 


[7] = 


I JO 

-+- 

p!i 


(II) 


{Y} = 


( 12 ) 


Applying the transformation (10) to the equation of motion (I), we obtain 
f7l T [MI ly] {Y} + [ 7 JT [C j [y] { Y} + [yf [K] fy] {Y} = [yjT {F} 


( 13 ) 


[Mr.jMr; 


1 • 

, X J | 

+ 

Cf.jCft- 
— — 

jx, | 

'XfiJKjV 
1 

x, j | F, + |3 t F 2 



xS] 


p* ip* 

V^2| '-22 

(xf| + 

_k?,!k 2 * 

XV j I f 2 


04) 


29 


where: 


fM,*,] = [M, ,] + j(J] T [M 2i ] + [M, 2 ] 131 +- £P} T [M 22 ] [p] 

= [M,,] + fp] T |M 22 ] 

IM*i) = |M 2 ,| + |M 22 J |p] 

[ M 22 1 = | M 22 1 

The partitions |C*| and [K*l have the same form as equations (15) through (18). That is 

|C,*,| = ie.il + ipf 1 |C,,| + 1C I2 | |pl + lp] T [C 22 ] |p] 

|K*| = |K,,| + lp| T |K :| | + |K I2 | (p| + |p| T [K 22 | |p] 

etc. 


05) 

(16) 

(17) 

(18) 


Since the {X.) displacements are known or prescribed, 
equation from (14) describing the {X 2 } variable, that is 


we need consider only the matrix 


|M 2 * 2 | {xj} + |C 2 * 2 | {XV} + |K 2 * 2 I {xS 


{F : } - [M 21 ] $,} - |C 2 *,] {X,} - [K*i] {X,} 


(19) 


We can expand equation (19), replacing the coefficient matrices by the relationships 
equations (15) through (18). Thus, equation (19) becomes 


as described in 


rM «l <*’> + ICj-’l f* 2 J + |K 22 | rxf} = IF,} - [M 2 ,l {X,( - t M 22 ) [f5] (X,) 

- IC 21 } {Xj} - (C 22 J [pi {X,} - (K 21 | {X,} - [K 22 ] [p] {x,} 

In view of equation (7) the term [K 22 | IP) {X,} appearing in equation (20) can be simplified as 


|K 22 ||p|{X,} « -|K 21 |{X,} . 


( 21 ) 


^matr“x.Te“ me <f " COnVeni,! " ce sake > ,ha ' ,he da "Pi»* ™«W*. (Cl. U pmportional to the .stiff- 


[C] = a [K] 


( 22 ) 


Substituting equat ons (21) and (22) into equation (20) yields 


[M 22 ] {Xj} + [C 22 ] {X c 2 } + [K 22 ] {X c 2 } 


{F 2 (t)} - [M 2 ,J {X,} - [M 22 ) [0] {X,} (23) 


nnnli P H T f h reSP f nS£ ° f th L d,s f retel y m °deled structure subject to base motion excitation, {X,}, and 
applied forcing function, {F 2 (t)}, can be described by solution of equation (23) and its auxiliary 

thetZ r ' S ,nt . erestin ^° l note the left side of equation (23) represents the structure with* 
he base motion coordinates {X,} constrained to zero. The base motion excitation {X,}, appears in 
the forcing function terms on the right side of equation (23). PP -' m 


Usually it is more convenient to solve the 
coordinates rather than the discrete coordinates of 
discrete coordinates and the normal coordinates is 


equations of motion in terms of modal or normal 
equation (23). The transformation relating the 


{Xj} = [<)>] { q } 


(24) 


where 


{q} — vector of normal coordinates 
I4>J = transformation matrix 


We choose |(|>) to be the modal matrix obtained from the eigenvalue analysis of 
IM22I {X 2 } + |K 22 ] {Xj} = {0} 

That is, each column of |(f>) is an eigenvector (mode shape) of equation (25). 

, + |T yie S | t ,,i ' U,in8 CqUa,i ° n (24) in '° eqUa,i, ’ n <23) and P^^uKiplying .he resulting equation by 
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[Meq] {q} + {2^j [Meql {q} 4- [<o 2 ] [Meq] {q} = (<|>1 T ({F 2 } - [M 2I J {X,} - [M 22 ] [0j {H,}) 

(26) 

where [Meql, {2£io], and [or] are diagonal matrices given as (Ref. Hurty & Rubenstein) 

(Meq] = (<|>] T [M 22 ] [4>] 

[2^o>] IMeql = (<|>] T [C 22 ] [<\>] 

(w 2 l [Meql_j= [cJ>] T [K 22 ] [<t>] 

. Now equation (26) is a very convenient form for numerical solution. In view of equations (8) and 
(24), the displacements {X 2 } can be determined as 

{X 2 > = iPl iXi) + [<j>] {q} (27) 

Thus, the response of the discretely modeled structure subject to base motion excitation can also be 
determined by evaluating equations (26) and (27). 
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